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We have studied spin-flip transitions between Zeeman sublevels in GaAs electron quantum dots. 
Several different mechanisms which originate from spin-orbit coupling are shown to be responsible 
for such processes. It is shown that spin-lattice relaxation for the electron localized in a quantum dot 
is much less effective than for the free electron. The spin-flip rates due to several other mechanisms 
not related to the spin-orbit interaction are also estimated. 
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I. INTRODUCTION. 

Quantum dots (QD's) in semiconductor heterostructures provide a unique opportunity to study the properties of 
the electron_jmiantum states in detail and manipulate the electrons in these "artificial atoms" in a controllable way 
(see reviewsoB) . The shape and size of quantum dots can be varied by changing the gate voltage. This also tunes the 
number of electrons in the dot. Besides, the electronic states can be significantly modified by a magnetic field applied 
C3 . perpendicular to the plane of the heterostructure. P 

Quantum dots are considered as possible candidates for building a quantum computer.!! The crucial point of the 
idea is the necessity to couple dots coherently and keep coherence on sufficiently long time scales. In this respect, 
there is a great demand in the theoretical estimation of the typical spin dephasing time of the electron in the QD. In 
our previous worko we have shown that the localized character of the electron wave functions in the QD's suppresses 
the most effective intrinsic spin-flip mechanisms related to the absence of inversion symmetry in GaAs-like crystals. 
This leads to an unusually low rate of spin-flip transitions. However, in Ref. ^ we concentrated on the case of inelastic 
transitions between the neighbouring quantized energy levels in the dot which corresponds to a relatively large energy 
transfer. On the other hand, the quantum bit was proposed to involve two Zeeman sublevels of the same orbital level. 
Therefore, in the present work we consider the transitions between such sublevels. Since the transition involves a 
t—{ ' fairly small energy transfer, the resuls are very different from those of Ref. ||. 

I/"") , As in Ref. |[ we concentrate on the spin-flip processes due to the spin-orbit interaction. This is the main source 
of the spin-flips for the three- and two-dimensional electron states in the GaAs-type crystal without an inversion 
center. Besides, in such a polar-type crystal one finds a strong coupling of electrons to the bosonic environment via 
the piezo-electric interaction with acoustic phonons. The combination of these two mechanisms provides an effective 
spin-lattice relaxation of free carriers in AjjjBy semiconductors and heterostructuresJj We show, however, that the 
. spin-lattice relaxation for the electron localized in the QD is much less effective. 

We have calculated the rates for the different spin-orbit related mechanisms which cause a spin-flip in the course 
of the phonon-assistant transition between the Zeeman sublevels. Besides, we have estimated the spin-flip rates due 
to several other mechanisms, for example, due to the fluctuating magnetic field produced by the fluctuating electron 
density in the leads or due to modulation of the hyperfine coupling with nuclei by lattice vibrations. 
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II. SPIN-ORBIT MECHANISMS. 



We consider the case of a strong confinement in the z-direction and the typical lateral dot size is of the order of 
thousand A and much larger than the width of the 2D layer in the z-direction. We begin with the following one-electron 
Hamiltonian that is derived from the Kane model (see Ref. ||) and describes the 2D electrons of the conduction band 
in the presence of magnetic field B, lateral confining potential U(r) and the phonons: 
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Here p = —ifr'V + (e/c)A is the 2D electron momentum operator, m the effective mass, <x the Pauli matrices. 
Axes x, y, z coincide with the main crystallographic ones with the z-axis along the normal to the 2D plane (the [100] 
orientation). The magnetic field has an arbitrary direction. The third term describes the spin- independent interaction 
with the phonons, including the piezoelectric ones. The fourth term is the Zeeman energy. The other three terms 
describe all possible spin-orbit effects. TLi stems from the absence of the inversion symmetry in the bulk.H Velocity 
j3 takes the values in the interval (1 -j- 3) • I0 5 cm/s for GaAs heterostructures. H.2 describes the spin-orbit splitting 
of the electron spectrum due to the strain field produced by the acoustic phonons. There <p x = {\/2){u X y,p y }+, 
<Py = -(l/2){ u yx,Px}+, <fz = (l/2){w za; ,Pa;jkf - (l/ 2 ){ M zy , p y }+ , where {,}+ denotes the anticommutator, u tj is the 
lattice strain tensor, and Vo = 8 x 10 7 cm/sE In GaAs the electron g-factor (g = —0.44) differs strongly from the free 
electron value go = 2 owing to the strong spin-orbit interaction which mixes the valence band and conduction band 
statcstl. The admixture depends on the lattice deformation. Coefficient g can be found within the Kane approach, 
g = {2m /V3m)(A/E !T l{d/E g ), d = -4.5 eV is one of the three deformation constants describing the strain effect on 
the hole band splittings, g » 10.4. 

The three terms in Hamiltonian ([!]) correspond to the three distinct mechanisms of the spin flip. The first mechanism 
is due to the spin-orbit admixture of state with an opposite spin. While without the spin-orbit interaction the Zeeman 
sublevels correspond to the orbital state with the spin up or down, the spin-orbit terms provide a small admixture of 
the state of the opposite spin to each sublevel. This enables the phonon-assistant transition between the two states. 
This mechanism corresponds to term Tix ■ 

The second and the third mechanisms are described by the H.2 and H.3 terms and correspond to two distinct kinds 
of direct spin-phonon coupling. Below we show that the admixture mechanism is actually a dominant one. 



A. Admixture mechanism. 



Let us show that for this mechanism the matrix element of the J7 p ^(r, t) operator for the spin-flip transition between 
the Zeeman levels is proportional to the product of the Zeeman energy and the phonon strain field. Since we deal 
with_a small energy transfer, we consider only the interaction with piezo-phonons, hence, for mode qa, a = l,t , we 
haveEl 



h/2ptu cia exp(iqr - i^ qa i)eA qa &q a + c.c.;A 



qa: 



qa 1 



(2) 



where the effective piezo-electric modulus Aq a of wave qa has been introduced, /3n~j is the piezo-tensor, £ = q/g, q 
the phonon wave vector, e the phonon unit polarization vector and p is the crystal mass density. For the crystal of 
cubic symmetry without an inversion center (class T^) tensor /3jjy has only non-zero components (all of them equal 
to each other) when all three indexes i, k,j are different, (3 xyz — j3 xzy = .... = /114. For GaAs ehu = 1.2 x l0 7 eV/ cm, 
see, for example, Ref. || 

The matrix element for the spin-flip transition between the Zeeman sublevels of orbital level n with emission of 
phonon qa is: 



ph 



E 



E n - E k - gp B B E n - E k + gp B B 



(3) 



where states n, k and corresponding energies E n , Ek are determined by first two terms in Hamiltonian |l]). The spin 
quantization axis coincides with the magnetic field vector. In the absence of a magnetic field the two terms in Eq.(|^) 
cancel each other since (Hi)^ = — (Wi)£„ and the matrix elements involving thp-.plionon operator are symmetric 
with respect to the interchange of indexes n and k. This "Van Vleck cancellation"ElJO is a consequence of Kramers' 
theorem and reduces the matrix element by a factor of gpsB /htoo, Hloq being the typical distance between the orbital 
levels in the dot. Note that this cancellation occurs for a spin-orbit Hapiiltonian of an arbitrary form. For instance, it 
could include the third order terms in the lateral momentum operator .□ This is in strong contrast with the cancellation 
of the linear in the j3 terms in the matrix elements for the spin-flip transition between different orbital levelsu, which 
results from the fact that spin-orbit terms TL\ are linear in the lateral momentum operators, p x ,y Expanding in the 
above formula with respect to the Zeeman energy, using relation (j>i) n k = (im/ ti) (E„ — Ek){xi) n k and the condition 
that the phonon wave length is much larger than the dot size (i.e. g^BB <C \/ms 2 HuJo, s is the sound velocity), we 
obtain the efective spin-flip Hamiltonian which acts on the subspace of the Zeeman sublevels of orbital level n: 
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where E XjV = —V Xyy U p h(x,y)/e is the phonon-induced electric field in the location of the dot. Here we introduce 
polarizability tensor a that may depend on B z . It is given by: 

( Tl ) I D \ n 2 (Xi) nm (Xk)mn (K s 

a ik ( B z) = -2e 2^ F _ F ( 5 ) 

Effective Hamiltonian (Q) is a very general one and can be used to calculate the spin-flip rates for arbitrary states 
and dots. To specify, we will consider only parabolic elliptic dots with the main axes along the x, y symmetry axes, 
u> Xt y being the oscillator frequencies. Then the symmetry of the kinetic coefficients ensures that (aiy + Oiyx ) = 0. 
We have to calculate the spin flip matrix elements < +1/2 | a x _ y | —1/2 > over functions = ±1/2, which are 

the eigenfunctions of operator u z > , where the z' axis is directed along the magnetic field vector. These functions 
are expressed through the eigenfunctions \ rn of a z operator: ^ = J2m=±i/2 4™ 2, *(f. ®i fyXm, where L^ 1 / 2 ) is the 
finite rotations matrbJil and tp, i? are the azimuthal and polar angles presenting B in the spherical coordinates. We 
substitute E x ^ y in terms of the boson creation/annihilation operators. Then for the square modulus of the spin flip 
matrix element that involves emission of a phonon with wave vector q we obtain: 
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+ a yy1y) g [("xxfe ~ a TO <2W cos 2^ - 2a xx a yy q x q y sm2(p\\ (6) 

The summing up over all q yields the rate due to the first mechanism 

2vr f d 3 q 



Tl = T I T2% £ c « I ^» (q«) I 2 " 

^ ^ ' a=Lt 



= 1 ggff (^r) [(aL + 0(l + cos^)-(aL-Osin^cos2^] + ^5) • (7) 

Here Ci — 1, C t = 2 and s/, s 4 are the longitudinal and transverse sound velocities. The anisotropy factors used are: 
^4q,z = 36/i 2 4 cos 2 sin 4 9 sin 2 cos 2 <j) where (f>, 6 are the azimuthal and polar angles of vector q. < A^ t >= 4ft, 2 4 < 
(6r£j/e z + £x£z e 2/ + Cyt^) 2 >= 2/i 2 4 [cos 2 sin 2 + sin 4 9(1 — 9 cos 2 9) sin 2 <f> cos 2 <^>] , where < ... > means averaging 
over the orientations of the e vector in the plane which is perpendicular to q. The averaging is done by the formula: 
< e^efe >= (l/2)(<5ifc — £;£fc)- As usual, in the case of finite temperature Eq.(^) should be multiplied by factors 
N u + 1(N U ) for the transition with emission (absorbsion) of a phonon, — l/(e huJ ^ T — l),fiu> = gfisB. Thus, in 
the case of high temperature T 3> g^sB the spin-flip rate will be proportional to (gfisB^T. 

In the particular case of a circular dot u> x = u y = u>q we have a xx (B z ) = a yy (B z ) = a xx (0) = e 2 /moj 2 . Then, for 
instance, for the transition between the Zeeman sublevels of the ground state of the circular dot with emission of a 
piezo- phonon we obtain: 

The dimensionless constant A p shows the strength of the effective spin-piezo-phonon coupling in the heterostructure 
and ranges from w 7 • 10 -3 to ~ 6 ■ 10~ 2 depending on (3. The spin-flip rate exhibits a very strong dependence on 
the Zeeman energy and lateral confinement energy loq. To give a number, Ti m 1.5 ■ 10 3 s _1 for Hljq — 10K and a 
relatively large magnetic field B = IT. 

Formula (|3|) is written with allowance for the wave function corrections of the first order with respect to the 
spin-orbit Hamiltonian. The corrections of the second order are described by the following spin-orbit Hamiltonian: 

Ti az = ^-S-z(xp y - yp x ) (9) 

Then, using this Hamiltonian in formula (|^) instead of "Hi, we can get a nonzero contribution to the spin flip matrix 
element even with zero Zeeman splitting in the denominator (but with taking into account the orbital magnetic field). 
Keeping again only the term which is linear in qr in the expansion of exponent exp(iqr), for the rate finally we obtain: 
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where 



D() = 2Re E x mm (L, + (eB,rV2c)) mn £> = _. % fd B\ 

E n -E m \ dy dx J 

In the absence of the magnetic field quantities D x , D y are identically equal to zero. Using the properties of the matrix 
elements for the linear oscillators we obtain that in the case of elliptic (circular) dots D x — D y — 0. Keeping the 
term which is quadratic in qr in the expansion of the exponent exp(iqr), we obtain a non-zero contribution but the 
corresponding spin-flip rate is smaller than contribution I\ by a factor of (/3/s) 2 (u> c /u>o) 2 < 1) here oj c = eB z /mc. It 
is also clear that in the case of irregular dots quantities D xy are not equal to zero. The ratio of the corresponding 
rate and Ti can be estimated as T 2 (mof3a/h) 2 , where r is a dimensionless parameter which describes the deviation 
from ellipticity and a is a dot size. Even when r ~ 1 this ratio is of the order of unity for a typical dot size a w 10 3 A 
Therefore, for t « 1 we can expect that contribution Eq.(]lO|) is much smaller than I\. 

Note that, besides term Tli which is linear in the 2D momentum, the initial Hamiltonian also contains the term 
which is cubic in the momentum: (l/2)a x {p x ,p 2 } + — (l/2)a y {p y ,p 2 }+. Again, in the presence of the orbital magnetic 
field we could get some contribution to the spin- flip rate. To this end, we need to calculate quantities D x , D y obtained 
from D x ,D y by replacing operator L z + (eB z r 2 /2c) by {\/2){p Xl p 2 } + or (l/2){p y7 p 2 } + . In the case of elliptic 
(circular) dots we obtain D x = D y = because of the symmetry. 

B. Direct spin-phonon coupling. 

Using the standard presentation for the strain tensor in terms of the acoustic phonon modes, we calculate the matrix 
element of H.2 for the electron spin-flip transition between the Zeeman sublevels of orbital state $ with emission of a 
phonon with momentum q: 

M T.i = x(2^J fee a + g a ej \ -{(p x + ip y ),exp(iqr)} + | $^ . (12) 

For simplicity, here we set B || z. Similar expressions were obtained in Ref. [K5] for a different problem. The total 
spin-flip rate is given by the Fermi golden rule: 

r 2 = x 1^ bB J ^(ql+ql)\^\l{(p x +^Py),eM^r)} + \^\ 2 5(h Sq -g^, B B). (13) 

The relevant phonon wave length is much larger than the dot size, which allows for further simplifications. We 
concentrate on a circular dot with confining frequency uj q . For the orbital states with n — and I = 0, ±1 (the ground 
and the first two excited states): 

r 2 = + T^fl 1 1 +i)!] 2 - (14) 

2407rps 7 ?T 2 v /wg + (o>2/4) 
The spin-flip rate produced by term H.3 does not depend on the structure of the orbital state and is given by 

.. (pj,g-B)'( Wg B) 3 . . 

rs = — ^v? — (15) 

Let us now compare the rates r 12 .3 obtained. All of them are proportional to the fifth power of energy splitting 
g/is-B, so that their ratio hardly depends on the magnetic field. First, we note that the ratio of T3 and T2 is of 
the order of (gs t /gVo) 2 « 7.8 • 1CP 3 <C 1. So that T 2 is more important. The ratio of Ti and T 2 is of the order of 
(ehi4/mVo/fi) 2 (rnf3 2 ms 2 /(hu}o) 4 )- For Hujo — 1 -7- 10K the ratio is of the order of 10 6 -j- 10 2 and increases only for 
larger dots that have smaller loq. Thus, we conclude that the admixture mechanism dominates. 

C. Two-phonon processes. 

The calculated rate Ti is small partly because of the small phonon density of the states at the scale of the Zeemaa. 
energy. On the other hand, for the case of the spin-flip transitions between the Zeeman levels of usual impuritytll 
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the two phonon processes under some conditions may become more important than the single phonon processes. At 
sufficiently small Zeeman splitting the contribution of the single phonon processes is very small, and with increasing 
temperature the role of the processes when one phonon is absorbed and the other is emitted is increased. It is also 
true for the case of a quantum dot and here we give some formulas which describe the contribution of such two 
phonon processes for GaAs quantum dots in several limiting cases. We also indicate the conditions under which these 
contributions can be important. 

If we treat the interaction with the phonons in the second order, we obtain processes in which a phonon is scattered 
from state p to state q while the electron spin flips. The effective matrix element contains transitions to an excited 
orbital state with the emission or absorption of a phonon and then transitions back to the ground state with the 
absorption or emission of a phonon. The spin may flip either in the first or second transition. The matrix element 
isllil: 

< V 2 >~ (^—\ (eh u ) 2 [N p (N q + If £ { [H ^ + f-* ] + [H -™ + ^1 1 , (16) 
\ps^/pqj \ -Aa-risq -A a + hsp 



H P, q = (*o",exp(ipr)*±)(*±,exp(iqr)* ), 

where N p is the Bose distribution function and A is the energy separation between the ground state whose wave 
function is Vl^o and the excited state whose wave function for spin up, say, is We can neglect the Zeeman energy 
in the denominators since no Van Vleck cancellation occurs here. We consider again interaction with piczo-phonons 
since deformation phonons become important at very high temperature (see below). For simplicity we consider here 
only the case when the magnetic field is perpendicular to the 2D plane and study the relaxation of S z spin component. 

As it was shown in Ref. |], there is a cancellation of the linear in the /3 terms in the matrix elements of type 
(\&o > ex P(*P r )^'a ) f° r the spin-flip transition between different orbital levels. This is a consequence of the fact 
that spin-orbit terms 7ii are linear in the lateral momentum operators, p x .y For that reason, quantities H pm are 
proportional to the first power of j3 only if one takes into account the Zeeman splitting in the electron spectrum.cl We 
consider the temperature in the interval gfJ-BB <C T <C hioo, where fkoo is the characteristic energy distance between 
the levels in the dot. Since hsp ~ hsq < T, then we can neglect the phonon energies in the denominators while 
calculating the contribution to < V2 > proportional to the first power of (3. It is apparent that the spin-flip rate has a 
different temperature dependence for the temperatures smaller and larger than To ~ ^/ms 2 hujo. At this characteristic 
temperature the phonon wave length is equal to lateral dot size A. For GaAs at hu>Q ~ 10K temperature To « IK. 
Let us give the estimate for the spin- flip rate in the case T « To, when pX,q\ <C 1. Here the estimate for H^„ is 
Hp q ~ (j3 1 \u>o)(gfj,BB /hujo)(Xq) 3 . Then the relaxation rate is: 

if ><r) . ,2./.) £ |< v 2 >? am* - ,) - msB] = fg ^y'' (§)' (17) 

In performing the integral over p, we have neglected gpsB in comparison to hsp. 

In the case T 3> To the momentum components parallel to the plane are estimated as qu A ~ 1, otherwise the matrix 
elements are exponentially small. As to the q z values, the contribution of the region q z ~ T/Tis 3> 511 is much smaller 
than that where q z ~ q» <C T/hs. Such is the case even without regard for the fact that for the [100] orientation 
of the 2D plane the effective piezo-modulus A qQ introduced above has a smallness q\\/q z <C 1. Thus, calculating the 
contribution from q z ~ gy <C T/hs and taking into account that N p = T /hsp S> 1, we obtain for the spin-flip rate in 
the case T < T < hu> a : 

T m (T)a ^i^Bn™'?>> (TV 



h P (fiu„) 7 /2 \T„) 

The contribution from the deformation phonons is much smaller. In the case T > To the characteristic q z ~ 
T/hs q\\ ~ 1/A, i.e. the deformation phonons are emitted almost perpendicular to the 2D plane. Then, for the 
deformation potential contribution to the spin-flip rate we obtain: ~ (A^/?i)(/3 2 /s 2 )(gpBB) 2 T 3 / (huio) 3 ms 2 , where 
Ad ~ (1 /2tt) (S 2 m 2 /ph 3 s) is the dimensionless constant which shows the strength of the electron interaction with 
deformation phonons. For GaAs A d w 10 -5 . Even at T ~ huj Q the value of T d /T^ ^ {hua / ms 2 ) 5 1 2 (A d f3 2 / k p s 2 ) 2 is 

much smaller than unity for any realistic hajQ. For example, at hujQ = 30K this ratio is ~ 0.03. 

(2) (2) 
Let us compare the two-phonon contribution with TiT/gfisB. We see that the two-phonon contribution rJ, 

prevails at sufficiently small Zeeman splittings: gpsB < ras 2 (A p s 2 / /3 2 ) 1 / 2 (T /To) 1 / 2 . Taking the maximal temperature 
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T ~ hujQi we obtain g/j B B < \(h p s 2 / I3 2 )yjms 2 /huvY^To. For ftuj « lOif we see that this contribution is more 
important for magnetic fields smaller than approximately 0.4 T (where the estimate for the spin-flip time is of the 
order of ms). On the other hand, at T ~ To we obtain for the same Hujq that g^sB < 0.03-K" (i.e. the two-phonon 
contribution is more important for magnetic fields smaller than w lkG). For these fields the characteristic spin-flip 
time is of the order of fs, i.e. it is still long. 

The general conclusion is that at sufficiently low temperatures (much smaller than hu>o) the characteristic Zeeman 
splittings below which the two-phonon contribution to the spin-flip rate dominates are small and corresponding spin- 
flip times are unusually long (see the estimates above). 

III. OTHER MECHANISMS OF THE SPIN FLIP. 

Let us discuss briefly other mechanisms of the spin flip. The spin transitions between the Zeepan sublevels of the 
impurity state in semiconductors ( mostly Si) were extensively studied quite a long time ago.l 14 rl 11 i Except spin-orbit 
coupling, several other mechanisms were proposed, such as: 1 ) modulation of the hyperfine coupling with nuclei by 
lattice vibrations, 2) the spin-spin interaction between the bound electron and the conduction electron in the leads, 
3) the spin-current interaction, when the bound electron spin flip is caused by the fluctuating magnetic field of the 
conduction electrons, 4) an exchange scattering process which flips the spins of both the conduction electron and the 
bound donor electron. Whereas the spin-orbit interaction strongly depends on the crystal symmetry and is. different 
for Si and GaAs, the other mechanisms are quite general in nature and we can profit from the discussion irjiSca. 

Mechanism (4) requires an overlap of the wave functions of the electrons in the leads and in the dot. In the context 
of QD it is considered in Ref. |l5|. The corresponding rates are not intrinsic to the dot since they are proportional to the 
barrier transparencies. They can be tuned to arbitrary low values. Refs.( [ll|Jl4] ) have demonstrated that the spin-flip 
rates associated with mechanisms (2,3) are very small. As an example, we give the rate estimation for mechanism (3). 
The Bio-Savar formula relates the magnetic field and current fluctuations in the leads so that < H 2 (l/c 2 a 2 ) < 

I 2 > u , a being the characteristic distance between the electron in the dot and the electrons in the leads. Using the 
Nyquist formula for the correlator of the currents we estimate < H 2 (l/c 2 a 2 )hu> coth(huj/2T)(l/R), R being the 

typical resistance of the leads or the dot environment. Thus, the corresponding spin-flip relaxation rate is estimated 
as 

r 4 ^ M |<ff 2 > w /fc 2 ^(^) 2 A (19) 

where A c = e 2 /moc 2 w 2.8 • 10 -13 cm is the classical electron radius, %u> = g/isB. Rate T4 is proportional to the 
first power of Zeeman splitting so that it may formally compete with Ti at sufficiently small splittings. However, 
this occurs at splittings that are so small that the corresponding rates are not observable. To give an example, we 
choose R = lOhm and Hujq = IK, which corresponds to a ~ 1.15 ■ 10 -5 cm. Then rate T4 would dominate if splitting 
hu> <C 2.5 • lO -3 !^. This corresponds to the rates lower than 8 • 10 _4 s _1 ! 

As to mechanism (1), i.e. modulation of the hyperfine coupling with nuclei by lattice vibrations, the relative strength 
of this mechanism and the spin-orbit interaction can be different for different materials. For example, in the case of 
Si where the spin-orbit interaction is much weaker than in GaAs, the dominant mechanism of the spin-flip for thp 
case of the phonon assisted transitions between the Zeeman levels of usual impurity (the situation studied in Ref.E 4 !) 
was found to be the modulation of the hyperfine coupling with nuclei by lattice vibrations. In the case of GaAs, 
however, our conclusion is that the dominant mechanism is the admixture mechanism of the spin-orbit interaction. 
This conclusion was reached for the first time in Ref. |l6|, where the calculations used essentially followed those in Ref. 
|l4] . Here we give the result obtained in Ref. [l6] for the rate due to the modulation of the hyperfine coupling with 
nuclei by lattice vibrations 

T h ~ {g^ B Bf 1 2 Lo 2 N /fi 2 s 5 p, (20) 

where ujn ~ (vqA 2 / a 2 zq) 1 / 2 /ft is the electron spin precession frequency in the random field of unpolarized nuclei, vq 
is the unit cell volume, A the hyperfine interaction constant, a the dot lateral size and zq the electron wave function 
extension in the z-direction. Finally, 7 ~ (l/m)(dm/dA) is the change in effective mass m with dilation, see also Ref. 
|l4| . For GaAs QD with a ps 10 3 A and zq ps 10 2 A, ojn can reach ~ 10 8 s _1 . Let us compare spin-flip rate Eq. (f20|) with 
Ti (this comparison was done earlier in Ref. |l^). Even taken for 7 ~ 50 (see alsca) we can easily see that Th will 
compete with Ti at the Zeeman splitting ~ ^uj^iTiLOQ) 2 / (3ehu which is so small that the corresponding rate is not 
observable. For example, for ftu>o = 10K the splitting is of the order of 10 -5 K. Therefore, the admixture mechanism 
of the spin-orbit interaction is the dominant one. 
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It should be noted, that in this work we have not considered the electron spin relaxation mechanism which is 
through the hyperfine interaction related to the internal nuclear dynamics.,— The latter is due to the dipole-dipole 
interaction between the nuclei which does not conserve the total nuclear spintZI. This mechanism might be important 
at low magnetic fields. However, this problem is not simple and needs a seperate investigation. 

Finally, we mention the experimental studies of spin relaxation in n-type GaAs quantum dots. Such an experiment 
has been recently carried outES. The non-equilibrium tunneling current through excited states in an AlGaAs/GaAs 
quantum dot was studied using a pulse-excitation technique which measures the energy relaxation time from the 
excited state to the ground state. Some excited states showed a relaxation time which was much longer than a few /is, 
while the other showed time much shorter than a few ns. This great difference in relaxation times was ascribed to the 
fact that some inelastic transitions are accompanied by the spin flip. For these transitions the relaxation time was so 
long that the method used in the above mentioned paper only allowed to give some estimation (much longer than a 
few /is). Though the transitions studied by T.Fujisawa et al. could in general involve the spin flip-transitions between 
the states with different orbital structures (this situation was considered in our previous papercl), the experimental 
data confirm the general statement that the spin-flip processes in n-type quantum dots can be really slow. 

In conclusion, we have calculated the rates for the phonon-assisted spin-flip transitions between the Zeeman sublevels 
in a quantum dot for all possible mechanisms and shown that the admixture mechanism of the spin-orbit interaction 
is a dominant one. The corresponding spin-flip rate Ti (see Eqs.(Q,^)) exhibits a strong dependence on Zeeman energy 
and at small magnetic fields takes very low values (up to seconds). 
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